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The most recent published report describing work funded by this grant is still Ref. 
the third such product based on the Doctoral Thesis of of Nitzsche (Ref. 2). Howevei, 
a paper by the Ph.D. candidate Valana Wells (Ref. 3) appears in the bound volume of 
the 14th I.C.A.S. Congress in Toulouse and was delivered with outstanding success by 
its author at that meeting. As described below for information, the Wells paper deals 
with a new linearized potential theory for a multi-bladed propeller in steady, subsonic 
flight. Travel expenses for this presentation were not furnished by NASA, nor was there 
an attribution to NASA because the research was supported form several sources. An 


extension is now planned, however, to the oscillating propeller for flutter analysis and the 
like; this will receive significant funding from the present grant. 

Professor Ashley described some of Nitzsche’s results in an invited seminar delivered 
at Arizona State University, Tempe, on March 30, 1984. Entitled “Vibration and Flutter 
of a Vertical Axis Wind Turbine Blade,” this not only reviewed Nitzsche’s theoretical 


investigation but also presented more general information about wind energy add about 



The September 1983 Progress Report gave some details about research undertaken last 
summer by the Principal Investigator and dealing with the influence on aeroelastic stability 
on various kinds of “strong nonlinearity.* During a busy academic year and summer, there 
has been time only for limited progress in this effort. Accordingly, a review of the subject is 
postponed until a later report. The remainder of the present one focuses again on the work 
of Ph.D. candidate John Green. Also described is a study of subsonic propeller theory by 
Ph.D. candidate Wells, whose Research Assistant appointment is being partially supported 
by Grant NGL-243. 

Activity of Graduate Student J. Green on Divergence of Laminated Composite Lifting 
Surfaces. 


The February 1984 Progress Report outlines the early stages of the work done on the 
divergence of forward swept wings. This early investigation focused on isotropic structures, 
for which analytic solutions were available (Refs. 4, 5 and 6). The integrating matrix 
technique proved both accurate and quick to converge upon the known results. 

The study of composite wings is well along, with comparisons being made to the work 
such as Blair, (Ref. 7). Blair’s paper includes both computational and experimental results 
for composite plates, with and without an aerodynamic shell, and it is the former case that 
is of most direct applicability. The largest part of Blair’s report deals with a composite 
plate with a fibreglass overlay to provide an aerodynamic shape. However, he also studies 
three cases with the plate alone and describes wind-tunnel tests for three fibre orientations 
and four sweep angles. Interestingly Ref. 7 concludes that tests done with the aerodynamic 
shell were unnecessarily complicated. 

The integrating matrix code developed at Stanford was run with these three fibre 
orientation angles and a full range of forward sweep angle. The calculations are summarized 
in Figs. 1, 2, 3, reproduced from the last progress report. The variable fibre angle 9 applies 
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to the two outermost fibres in the following lay-up: 


[T45/90/ =F 45/90/ 45/9O 2 /0 2 ], 

The figures show that the integrating matrix results match the experimental ones very 
well and that they are consistently closer than those from the OWING code which Blair 
used, This gives encouragement for the potential usefulness of integrating matrices as a 
fairly simple and cheap solution technique, since the data in these figures derive from the 
solution only of a 4 x 4 eigenvalue problem. 

The remaining Figs. 4 and 5 demonstrate two methods of displaying the behavior of 
composite wings with varying fibre orientations to a designer. The contour plot is more 
familiar, and permits more precise readings to be taken, but the three dimensional plot 
provides a very graphic demonstration, and is better for giving a feel for the phenomenon. 
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Fig. 4 — 3D Plot of Divergence Dynamic Pressure vs wing sweep angle and fiber orientation. 
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Before going on to complete the study of divergence as a function of the fcregoing 
and other parameters, it is useful to be able to evaluate the effect of flutter. Although 
divergence has been seen as the primary barrier to practical application of forward swept 
wings, the cure afforded by aeroelastic tailoring comes with a penalty. The composite 
layups that alleviate divergence can also reduce the flutter speed. (See for instance Refs. 
8, 10 and 11.) 

The aim here is to be able to modify the 3-D plot for divergence (e.g., Fig. 4) to 
include flutter and thus produce a general aeroelastic stability boundary. The method 
used to calculate the flutter speeds is based, of course, on the integrating matrix, and an 
outline of the relevant theory follows (Ref. 9). 

The general vector equation of motion is 

£ = Zy + s 1 My + sCy-Q(,,,)y, (1) 

where 

y = state vector 
Z — structural matrix 
M = mass matrix 
C — damping matrix 
Q — aerodynamic matrix 
s = Laplace variable 
q = dynamic pressure 

After applying the integrating matrix to both sides of the equation (1) and enforcing a 
cantilever boundary condition, the problem takes on the following matrix form: 

( I + Gli) G 2 S C?24 

G32 (I + (733) f ?34 

Gi 2 O43 (/ -f O44) 

with 

G = T [Afa* 2 — <?(s‘,A)] (3) 
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Here 

T =■ includes boundary conditions, structural matrix and integrating matrix 

3 * = non dimensional Laplace variable 

A=nondimensional dynamic pressure 

Equation (2) depends upon the two variables, and flutter is indicated when the real 
part of s* becomes positive. Thus to solve (2) it is necessary to assign a value to A and 
calculate the value of s * for which the determinant is zero. This solution is done using 
an algorithm similar to Muller’s method (Ref. 9) to find the roots of the determinant. 
Muller’s method proves to be a very efficient and accurate approach. 

Having found the root of the determinant, a check is made to see if its real part is 
positive. If not, a new value of A is tested. When the real part becomes positive a further 
calculation can be made to evaluate the value of A at which flutter occurs. 

In order to be able to generate a 3-D carpet plot like Fig. 4 it is necessary to calcu- 
late the flutter speed of a large number of configurations. If the flutter speed has to be 
determined by manually testing the sign of the root then obviously much time would be 
wasted. Thus it is planned to develop an automatic “search” routine such as was used by 
Edwards (Ref. 12). He applied a gradient search algorithm with some success, and it is 
hoped that such a scheme will make 3-D plots feasible. 

The results of the divergence calculations were compared with the experimental results 
of Blair (Ref. 7); however, he did not include any flutter examples. It is intended to 
compare the flutter calculations with the experimental work of Sherrer et.al (Ref. 8), 
but since all their tests were for a composite plate with an aerodynamic shell, it will be 
necessary to include a correction for such a shell in the integrating matrix calculations. 
This work is currently being undertaken. 


Activity of Graduate Student V. Wells on Subsonic Propeller Theory, with a Goal 
of Aeroelastic Analyse 

This summary of progress is baaed primarily on the theory and examples of Ref. 3, 
copies of which can be supplied if required. After reviewing past linearized solutions for 
incompressible and compressible flow past thin screw propellers, the wo7k focuses on a 
potential flow solution in the reference frame depicted by Fig. 1 of the February 1984 
Progress Report. Writing a linearized P.D.E. for the disturbance potential, the procedure 
converts to a non-inertial cylindrical system (r,0,z,t) defined by 


r = r, 9 = 0 + ut 
z — z, t = t. 


( 4 ) 


Here oj is angular speed in md./sec. Performing these coordinate transformations gives 
the unsteady, linearized potential equation in the rotating system: 
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( 5 ) 


where ia the radially varying apparent Mach number in the angular direction. 

In general, the solution to a linearized lifting body problem is subject to the following 
conditions: 

1. ) Tangency of the flow at the body surface; 

2. ) Vanishing (or at least finiteness) of disturbances an infinite distance away from the 

body and its wake; and 

3. ) Zero normal velocity through the wake. 

The first of these conditions has the mathematical representation 


dF - 

~ + V • VF = 0 on F = 0 
dt 




( 6 ) 


tZi: J* 
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where F = 0 represents the surface of the lifting body, and \? = V# is the total velocity 
vector in an inertial coordinate system. By the rule for vector time derivatives, in a non- 
inertial frame, the total velocity, seen in the coordinate system rotating with angular 
speed, (j, is: 

W = V ~tixr=:P~ur0. (7) 

Substitution into the above boundary condition, along with the conversion of the time 
derivative, giveti the boundary condition hi the rotating reference frame: 

°"W~‘‘’W + ^'' VF = 17 + ^' Vf «nF = 0, (8) 

which is the expected result. In other words, the body boundary condition may be treated 

in the same manner in the rotating coordinate system as it is in the inertial system, as 

*** 

long as the local freestream velocity, W, is considered. 

The wake boundary condition can be made mathematically tractable by assuming - 

1. ) Unsteady perturbations are small so that the unsteady wake has, to leading order, 

the same shape as the steady wake; and 

2. ) The slipstream velocity, w, is much smaller that the freestream velocity, U. 

With these assumptions, the position of the wake can be approximated by: 

9 — ~z — +(n — 1)—-, r<R> ’ (9) 

for B = the number of propeller blades, n s the blade index, and R = the blade radius. 
For r = ro and n = 1, for example, the equation describes a single helix with radius r 0 
and origin at 9 = 0. The wake surface is a helicoid made up of an infinite number of these 
helices with radii from r = 0 to r = R originating at z = 0 and at each angular station 
corresponding to a propeller blade location. 
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The first order analysis assumes that all motion of the helicoidal wake occurs in the ax- 
ial direction. Since there can be no flow through the wake, the wake condition is expressed 
as 


UJ COS € = u t COS € - Sill € 


( 10 ) 


on 

= r < R (11) 

where u, s induced velocity in the z-direction, u 6 = induced velocity in the 0-direction, 
and e = helix angle = tanr'U/ur. In terms of the velocity potential, then, the boundary 
condition is 


dip 

txiroj ss ujr~~ - 
dz 


Udp 

r 30* 


(12) 


Most work to date has focused on the aerodynamics of the steady propeller; that is, a 
rigid propeller experiencing a uniform inflow. These conditions allow the time dependent 
terms in equation (5) to be excluded since, now, the flow field appears steady in the 
rotating, translating coordinate system. 

The dimensional spatial variables, r and z, have convenient dimensionless forms — 
p — ur/U and 5 = Dropping the bar over the dimensionless z variable, the steady, 
potential equation becomes: 


dp! pdp P 


dz* 


+ -i(l-MV) 


d 2 <p 

~df 


2 M 2 


d 2 tp 

d&dz 


= 0 , 


(13) 


with (3 2 = 1 — M 2 . Equation (13) appears to depend on one parameter (M) only. However, 
the advance ratio parameter is hidden in the dimensionless variables, p and z. 
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In “helical" coordinates described by £ = 6 - z and cr = 0 + *, this becomes 

W + -JP* + (1 + + (1 + ^ - 4M*)p<„r - 2(1 + ~ - 2M*)<p (ff = 0. (14) 

Goldstein (Ref. 15) shows that, in the far wake, the flow field depends only on the variables 
p and £ . (^ runs counter to a helix at a given radius.) This results from purely geometrical 
considerations, and therefore holds for both compressible and incompressible flows. Thus 
in thfi far wake, 

<Pp, + -<p, + (1 + ~ W = 0, (15) 

which is exactly the equation solved by Goldstein. This result is of interest on two accounts: 

1) The perturbed flow in the far wake has an incompressible character which implies that 
quantities which depend only on the far wake are invariant from incompressible to 
compressible flow. 

2) The far wake solution for a propeller in incompressible flow also satisfies the compress- 
ible equation. 

Davidson (Ref. 16) pointed out the latter and used that fact in his study of the 
propeller in compressible flow in a wind tunnel. This study also utilizes the far wake 
solution as will be shown. 

The problem of compressible, potential flow about a propeller reduces to that of solving 
a partial differential equation Eq. (14), subject to 

1. tp finite at p = 0, 

2. Sommerfield radiation as p becomes infinite, 

3. No flow through the propeller blades, and 

4. No flow through the propeller wake. 

The solution to the equation for the incompressible velocity potential, <pi[p,$-z), also 
solves the compressible equation. Though it cannot capture the compressible effects at the 
propeller blades, this solution accurately describes the far wake flow. Thus, a superposition 
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of a compressible solution, which accounts for effects near the blades but dies out ii% the 
far wake, on the incompressible solution also constitutes a solution to the linear differential 
equation. As will be shown later, the superposed solution is desirable in, this case as the 
analytic form of the compressible solution does not lend itself well to satisfying the wake 
boundary condition. 

Though most forms of the incompressible velocity potential, <pj f are derived at an 
infinite distance behind the propeller plane, any form given as <pi(p } 9 - z) must solve 
Eq. (14). However, in order for the superposed solution <p = <p e 4- <p> to form a solution to 
the entire problem, the sum, <pi + <p ei must satisfy the boundary conditions. The first two 
conditions, finite <p at p = 0 and proper dying out of the disturbance as p — ► oo, do not 
cause any difficulty as they can be satisfied independently by ipt and <p e . 

In order to account for the more problematic conditions (3) and (4) above, the flow field 
can be artificially divided into two regions - one forward of the propeller plane in which only 
the compressible velocity potential holds, and the other behind the propeller plane in which 
both <pi and <p e apply. The form of <pi , which automatically provides for the existence of a 
trailing wake, necessitates the division of the flow field. This division, however, requires the 
introduction of two more “boundary” conditions which assure continuity at the propeller 
plane. Designating <p e + as the forward potential and <p~ = ( <pi + ip e ~) as that behind the 
propeller, and situating the propeller at z = 0 as before, a continuous flow field is ensured 
by 


<Pc 


+ 





(16) 


and 


d<p e + j 

_ 9(pi j 

9(p e ~ 

dz ' 

s=o dz '* = o 

dz 


(17) 


Since both <p + and <p~ are valid at z — , 0, both must satisfy the condition of no flow 
through the propeller blades. 
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Details of this incompressible result are worked out in Ref. 3 on the basis of Reissner’s 
paper (Ref. 17). A general solution to the steady, compressible potential equation may be 
obtained through separation of variables. The resulting series solution has the form: 

v = £ A n u , p ) U ^" B, \ (18) 


Up) = c-fl (esjMPrfB* + a;) , 


(19) 


and 

n 2 B 2 M 2 ) 

pi /*' 

C n b represents a linear combination of solutions to Bessels equation, A„ is the separation 
constant, or eigenvalue, and p, ss f3jM is the dimensionless radius at which the local 
freestream velocity becomes son'ic. Applying the boundary condition requiring finiteness 
of ip at p = 0 gives 



Z n {z) =5 exp 


.nB , 1 

^ ± - 
ps 2 




■^n(p) = JnB (psJW-vPB* + Xlj . (21) 

A separated solution is a slightly different form from that given above was first proposed 
by Busemann (Ref. 13). 

The behavior of the J-Bessel function as its argument approaches infinity will cause 
automatic satisfaction of the radiation boundary condition. This leaves no mechanism for 
determining values for the eigenvalues, A„; in fact, the notion of the eigenvalue has no real 
meaning for problems of infinite domain. Therefore, the solution must be obtained using 
a method other than a simple separation of variables. 

Due to periodicity in the angular coordinate, the expected solution to the partial 
differential equation has the form: 

<P=fl ( Pn(p,zy nB9 ( 22 ) 

n=0 
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Substituting this into the PDE gives 


0 ■ 
05 - r >/-i : 


+ j<Pn, ~ (l + ^r» = -p 2 <Pn„ + 2i»SAf*^ Bf , (23) 

Performing a Hankel transform on this equation gives 


f 3 2 ip„ - 2\nBM*V t - ( 7 * - i» 2 £*Af*) p = 0 , (24) 


a second order, ordinary differential equation for the Hankel-transformed potential, ip. The 
Hankel parameter, 7 , appears only parametrically. The solution to this equation is easily 
obtained: 


<P„(z\ 7) = A n exp 


{ 


.nB . 1 

_ n 2 B 2 l 

[ , 7 i ± n 

w • 

1 



Inverting this transformed solution results in 


(25) 



P )cxp 


.nB , 1 

. n 2 fl 2 

t — — ± — . 

p\ 

1 

w » 

l 

•4 

t- 



dl- 


The complete solution to the compressible potential equation is then 


(26) 


n=0 *'° 


An(7)7^r»B('7P)e®P 


.nB , 1 

2 n 2 B 2 

'IT P\ 

l 2 2 

P 2 , 


d7. 


(27) 


It remains to determine the unknown “constants,” A„( 7 ). 

A careful “matching” process described in Ref. 3 makes it possible to relate all of the 
required constants to the single unknown function V(p), which is the unknown distribution 
of circulation along each of the identical propeller blades. Table 1 summarizes the results. 
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Table 1 — Coefficients for Compressible Velocity Potential Series 

After the convergence of integrals appearing in Eqs. (26) and (27) is proved, one 
proceeds to apply a flow-tangency boundary condition to determine I\ This condition 
may be expressed as 


u 2 cose- ursine — Wa for z,0 -» blade, (28) 

for the simple case of uncambered, thin airfoils sections. u t and u* are the perturbation 
velocities, achieved by differentiating the velocity potential. u,,U0,£,W, and a are all 
functions of radial position. 

In the lifting line approximation, which is developed in detail, there is no blade at 
which to apply condition (28). Instead, the condition given by Reissner must be used. 
Referring to figure 6, this condition is given as 

a = a a - ~ at z 0, 9 (29) 

a g is the geometric angle of attack, and t /»• is the induced velocity perpendicular to W. 
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Figure 6 — Velocities at the Blade 


Two well-known relations assist in reducing this to usable form: i — pW T and i — 
%W 2 cct a a, where l = lift per unit radius, n = density, q a = section lift curve slope and 
c = blade chord. Combining these two expression gives a = 2Y jW cc\ a . Substituting this 
into (29) gives 


Vi 2VU 
W + p a cc, a W ~ a °' 


z — ► 0, 9 — ► 


7T 

I’ 


(30) 


where c = c/R , an inverse “blade aspect ration.” Note that — u : cose — uesine, as 
before. 

Condition (30) introduces two parameters, c and Cf a . c is easily estimated for a typical 
propeller blade; however, ci a is not known in advance for compressible flow. It is assumed 
here that an estimate of c/ a based on strip theory will introduce only minimal error. 

By defining <p* = ^<p and T* — ftT, Eq. (30) can be written in dimensionless form: 


d<p* 

p ~d. 7 


1 d<p* 
p d9 



= {l + p 2 )a g ,z 


0 , 0 — ► 


B' 


(31) 
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Carrying out and collecting the expressions for the derivatives in Eq. (1)1) leads to a 
rather complicated integral equation pair for the circulation distribution (see Eqs. (4.5 
a,b) of Ref. 3). A useful method for solving these relations involves choosing a series 
“shape function” for T*(p) with unknown coefficients. Collocation at a set of K stations 
along the blade then yields K X K linear equations in K unknown coefficients. 

The first attempt at selecting these shape functions followed a suggestion of Reissner 
(Ref. 17). This led, however, to difficulties with the uniformity of convergence and has been 
improved upon since Ref. 3 w^s written. At Mach Number 0.65, for example, satisfactory 
convergence is now attained with approximately K = 10 on a typical two-bladed propeller. 
It turns out, incidentally, that no simple “compressibility correction” can be proposed by 
analogy with the Prandtl-Glauert correction for steady flow over subsonic wings. 

Figure 7 compares the compressible and incompressible circulation distributions for 
propellers with the same geometric angle of attack distribution. As expected, the com- 
pressible case has a greater magnitude of circulation, paralleling the same trend due to 
Mach number which occurs for straight wings. 



p/o« 

Fig 7 — Circulation Distributions 
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Favorable comparisons have also been made with test results on two-bladed machines, 
when an empirical adjustment is made to account for the profile drag of blade sections. 

The results obtained inspire confidence in the derivation and the solution methods 
since a) they show convergence in both the series solutions for the velocities and the shape 
function solutions, and b) they show trends anticipated from wing theory. Of course, a 
favorable comparison with experiment would reinforce confidence in the method. 

Further examples can be readily treated and will be shortly forthcoming. These include: 

a) thrust and torque distributions, 

b) efficiency calculations, 

c) effects of changing Mach number and advance ratio, and 

d) unsteady motions. 

The Weissinger boundary condition (Ref. 18) will also be studied for accuracy improve- 
ment. In addition, a viscous drag correction term must be included to correctly predict 
the thrust and torque distributions, particularly for a compressible propeller analysis. 

Activity of Graduate Student In Lee on Cross-Sectional Resonances of Wind Tunnels. 

A limited description of this research was provided in the June 1984 continuation 
proposal (Ref. 18) for the present grant. Extensive details are not repeated herein. Mr. 
Lee’s work has important implications for high subsonic flutter testing, but it has received 
in the past only a small amount of support from the grant. Whether it will be possible to 
continue this funding in the future depends, of course, upon other demands on the limited 
available funds and on completion dates of other student research. 

It has been well known since two-dimensional analyses at NACA Langley in the 1950’s 
and subsequent Brittish studies on 3-D cross-sections that there are various potentially- 
significant wall effects on the ^irloading of wind-tunnel models of lifting surfaces performing 
unsteady motions in a subsonic airstream. ( In the interests of brevity, the many relevant 
references are not cited here.) Although transonic effects can be significant, probably 
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the most dramatic of these occurs when a frequency of oscillation “resonates” with a 
natural acoustic frequency of the tunnel cross section itself. Such “resonances” are most 
pronounced in closed tunnels and when the fundamental tunnel mode is involved. Effects 
on flutter speed of the order of 20% have been demonstrated experimentally. Thus it is 
desirable to be able to predict these modes and frequencies accurately. 

Although the theory can be developed in a more sophisticated fashion, the key idea 
is that the “resonances” consist of 2-D standing waves in the cross section, with an effec- 
tive speed of sound aoo\A — M 2 . Here a is the undisturbed free-stream acoustic speed 
and M < 1 the test Mach No. It is significant that this effective speed, and the associ- 
ated frequencies, become quite low as the transonic range is approached. (Ultimately, of 
course, nonlinearity and the formation of shocks tend to invalidate the linear potential-flow 
idealization.) 

What Mr. Lee has succeeded in doing is to program an accurate, convergent finite- 
element code for finding 2-D acoustic modes of an essentially arbitrary tunnel shape. 
He has analyzed, for instance, the octagonal configuration of the NASA Langley TDT, 
discovering somewhat different natural frequencies from those previously obtained from a 
rectangular approximation to its shape. His program, in the original version, was capable 
of handling open or solid wall boundaries, along with combinations of open and solid 
portions that could be accomodated by the finite-element formulation. Many examples 
have been calculated. 

Mr. Lee is currently engaged in examining porous and slotted boundaries of the sort 
employed in transonic facilities. Following up on some British work, he has also discovered 
interesting influences due to a finite plenum chamber (at constant pressure) surrounding 
such a transonic test section. Publications are anticipated in the near future. 
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